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Abstract 

We present a unified treatment of the phase space of a spatially flat homo- 
geneous and isotropic universe dominated by a phantom field. Results on the 
dynamics and the late time attractors (Big Rip, de Sitter, etc.) are derived with- 
out specifying the form of the phantom potential, using only general assumptions 
on its shape. Many results found in the literature are quickly recovered and pre- 
dictions are made for new scenarios. 
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1 Introduction 



There is now substantial agreement that the expansion of the universe is accelerated, 
with evidence from type la supernovae [I], WMAP data |2j, and large scale structure 
surveys Two classes of models aim at explaining the observed cosmic acceleration: 
one modifies gravity on large scales by introducing corrections to the Einstein-Hilbert 
[3] or the Palatini [S] action (often there are instabilities jH] or incorrect post-Newtonian 
limits [7]), and the other advocates the existence of dark energy with density p and 
exotic pressure P < —p/3. 

Most dark energy candidates proposed thus far are scalar fields; there is marginal 
evidence |J] for an effective equation of state parameter of the dark energy w = P '/ p < 
— 1, which is equivalent to increasing Hubble parameter H > 0. If confirmed, this 
measurement is important because an increasing H cannot be explained by Einstein 
gravity with a canonical scalar field [B1E1E]- If the universe really superaccelerates 
it may end its existence in a finite time in an explosive expansion of the scale factor 
accompanied by diverging dark energy density, called a Big Rip [T21 E| • Other kinds of 
singularities at a finite time in the future discussed in the literature are called "sudden 
future singularities" and are characterized by finite scale factor and Hubble parameter, 
and diverging H and dark pressure [Tlj . 

Simple models of a superaccelerated universe employ a scalar field coupled nonmin- 
imally to the Ricci curvature PHI El EH El E] or a phantom field, i.e., a minimally 
coupled scalar field with negative kinetic energy (HI El E] ■ Phantom fields and fields 
with non-canonical energy are present in string theories jTHl El 120] and supergravity 
[2*T] . are associated to bulk viscosity due to particle production [22], or arise in higher 
order theories [23 . It has also been proposed that early inflation and late time accel- 
eration of the universe can be unified in a single theory based on a phantom field [23] , 
A phantom field poses several challenges: it is subject to severe instabilities, which may 
perhaps be avoided by thinking of the phantom as an effective field theory resulting from 
some fundamental theory with positive energies [21]. A fundamental quantum phantom 
is very difficult to stabilize [T31 E] ■ 

In general, a form of dark energy with w < — 1 is worrisome because it violates 
the energy conditions cherished by most physicists and opens the door to disturbing 
exotica. For example, even a small amount of exotic matter violating the weak energy 
condition may be sufficient to open up a wormhole and make time travel possible |27j . 
On the other hand, even a simple classical scalar field coupled nonminimally to the Ricci 
curvature may violate all of the energy conditions [2H] , and one should probably not be 
too rigid in rejecting phantom fields a priori. 

The peculiar dynamics arising from a negative kinetic energy density has been stud- 
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ied with the help of a toy model consisting of two mutually coupled oscillators, one 
with negative kinetic energy representing the phantom field, and the other with pos- 
itive kinetic energy representing the gravitational field [T^]. However, this toy model 
fails to properly describe the dynamics [29, and it is preferable to study the actual 
dynamical system. Furthermore, a number of papers in the literature study phantom 
models with different specific potentials [SO]- Here we study the dynamics without toy 
models or approximations and we deduce the behaviour and the late time state of a 
phantom-dominated universe without specifying the form of the potential. The assump- 
tion of a spatially flat Friedmann-Lemaitre-Robertson-Walker (FLRW) universe plus 
general assumptions on the potential (such as boundedness, monotonicity, etc.) allow 
one to derive the asymptotic dynamics. In addition to the need to understand phantom 
dynamics with general potentials, our study is motivated by the question whether the 
universe will end in a Big Rip or will expand forever. In fact, a universe that super- 
accelerates may end its existence in a Big Rip or it may expand forever, depending on 
the dynamical equations. Often authors working with a specific choice of the phantom 
potential find a de Sitter regime as the final stage of evolution (late time attractor). We 
show that certain features of the phantom potential determine whether the final state is 
an indefinite expansion or a Big Rip, and the analysis needs not be repeated for all the 
possible potentials. We adopt the notations and conventions of |32j and units such that 
c = G=h=l. 

2 The phase space of phantom cosmology 

Based on the recent cosmological observations we adopt the FLRW line element 

ds 2 = -dt 2 + a 2 {t) (dx 2 + dy 2 + dz 2 ) (2.1) 

in comoving coordinates (t, x, y, z). We consider the situation in which dark energy has 
already come to dominate the cosmic dynamics, hence a phantom field is the only form 
of matter in the Einstein equations. The energy density and pressure of the phantom 
are, respectively, 

P=-^0 2 + V(0), P = -i0 2 -V(0). (2.2) 

The distinguishing negative kinetic energy is evident in eqs. ()2.2|) which exhibit the 
"wrong" sign for the kinetic term 4> 2 /2. The Einstein equations for and the scale 
factor a{t) are 

H 2 = 1 (-^ + 2V) , (2.3) 
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H + H 2 = - (0 2 + V 
dV 

+ 3H<p - — = , 



(2.4) 
(2.5) 



where H = a/a is the Hubble parameter, an overdot denotes differentiation with respect 
to the comoving time t, and k = 8nG. The sign of the term dV/d<fi in the Klein-Gordon 
equation ()2.5|) is the opposite of the usual one. The potential is required to be positive 
by eq. ()2.3j) . and this modification of the usual Klein-Gordon equation implies that a 
phantom falls up in an increasing potential, is repelled by a minimum, and settles in a 
maximum of V |12M33| l3*4]. In addition, a monotonically increasing unbounded potential 
may be expected to generate runaway solutions - the analogy with a canonical scalar 
field in a negative potential not bounded from below applies (these properties are proven 
later in this paper). 

Eq. (J2.3)) implies that the phantom energy density is non-negative for any solution 
of eqs. (|2.3j) - (|2.5j) . Only two equations in the set (|2.3|) - (|2.5j) are independent - when 
7^ one can derive the Klein-Gordon equation ()2.5|) from the other two. 

The field equations can be obtained from the Lagrangian L or the Hamiltonian 7i 



L = 3aa + na 



+ V(<f>) 



n 



H 2 + 



K 

6 



kV (</>)' 



(2.6) 



the Hamiltonian constraint (J2.3|) corresponding to 7i = 0. The orbits of the solutions 
of eqs. (|2.3|) - (|2.5|) are constrained to the surface of constant energy Ti = in the phase 

space ^a, a, 0, tpj . We choose as dynamical variables the Hubble parameter and the 

scalar field (H, (f>) which correspond to physical observables 1 . The phase space accessible 
to the orbits of the solutions is the two-dimensional energy surface 7i = in the three- 
dimensional space ^H, 0, (fij ; this surface is in general curved (the appendix provides an 
example). 

Once the values of H and are given, one computes the corresponding value (s) of 
by rewriting the Hamiltonian constraint ()2.3|) as the algebraic equation for 



2V + 



6H 2 







/v 



(2.7) 



1 Other authors choose as dynamical variables various combinations of H and (f> but this obscures 
the physical interpretation of their results. 
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with solutions 

3H 2 




v{4>) 



K 



= ±y/F(H,<l>). (2.i 



Depending the form of V ((f)), there may be a region T of the phase space forbidden to 
the dynamics and corresponding to a negative sign of F = 2 (V — 3H 2 /k): 

T=^{HAA)-- V(0)<^}. (2.9) 

When F (H, 0) > there are two distinct values of for each value of H and 0, corre- 
sponding to the fact that the two-dimensional phase space consists of two sheets joining 
each other only at the points of the set 

B = {(H,<P,£): K(0) = ^,0 = oJ, (2.10) 

which constitutes the boundary of the forbidden region F and lies in the (H, 0) plane. 
We denote the sheet corresponding to the positive sign of "upper sheet", while the 
other is called "lower sheet" . This structure of the phase space is general in scalar-tensor 
cosmology 31 and is well known in the theory of a scalar field with canonical kinetic 
energy coupled non-minimally to the Ricci curvature (UHl EEl EZ] . Phantom cosmology 
can be seen as a special case of scalar-tensor theory for which there is always a forbidden 
region in the phase space, while this needs not be the case in other scalar-tensor theories. 
The equilibrium points of the system ()2.4|) - (j2.5|) consist of de Sitter spaces with constant 
scalar field (H Q , O ). If they exist, these fixed points satisfy the constraints imposed by 
eqs. (I2~3|) - (123J) 

H = ±^, K = 0, (2.11) 

where Vq = V((f)o) and Vq = ^jr . The Hamiltonian constraint (|2.3J) can only be satis- 

fied if V(0) > 0, which guarantees that H is real - then there are de Sitter fixed points 
provided that dV/ d<p has zeros (the existence of equilibrium points is not guaranteed in 
all scalar-tensor theories). These fixed points lie on the boundary B of the forbidden 
region 2 . 



2 In other scalar-tensor theories the fixed points, if they exist, may be located anywhere on the energy 
surface H = 0. 



5 



What happens to an orbit with initial conditions chosen exactly on B ? To answer 
this question we consider the tangent to the orbits 

f = (H, 0, 0) = (J 2 , 0, ^ - aff^ (2.12) 

in the ^if, 0, 0^ space. On $ it is = and T& = (0, 0, dV/d(p). Hence if dV/d(p\ B > 0, 

an orbit beginning on B will move into the upper sheet, while if dV/dcj)\ B < 0, the orbit 
will move into the lower sheet. If instead dV/d<p\ B = one has a fixed point satisfying 
the properties H = 0, = 0, H 2 = kV/3, and V = 0. 

The following result is true for any form of the potential V((p): H always increases 
along the orbit of any solution except at points where = 0, at which also H vanishes. 

This follows by combining eqs. ()2.3|) and ()2.4j) to obtain 3 

H = ^\ (2.13) 

hence H > except at the points where = 0, in particular the fixed points (H , O ). 

Eq. (I2.13J) is expressed by saying that the universe always superaccelerates (i.e., H > 
as opposed to acceleration defined by d = H + H 2 > 0), or that the phantom field is a 
form of superquintessence [TTH ITT] . Eqs. (J2.3)) and (J2.13j) imply that 

H + 3H 2 = kV . (2.14) 

A second result is that there are no limit cycles (periodic orbits). 

In fact, H(t) is non-decreasing and hence it cannot periodically return to its initial 
value apart from the trivial case of a fixed point. If dV/dcf) has definite sign this possibility 
is also excluded, cf. eq. (|2.11j) . 



3 Bounded potential 

If the potential V(<j)) is bounded from above by a (positive) constant Vq, then the asymp- 
totic solution of eqs. at large times is such that H — > and — > 0. The 
resulting de Sitter attractor (H , O ) is a global attractor. 

In fact, according to eq. ()2.14|) . H = kV — 3H 2 ; assume that H does not tend to zero 
(or that it does not tend to zero faster than t -1 ) as t —>■ +00. Then the limit 

lim H{t) = lim / H{t) dr = +00 (3.1) 



5 Eq. (|2.1^(l is equivalent to eq. (7) of ^TJ multiplied by 
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since H > 0. But then V(t) > 3H 2 /n — > +00 as t — > +00 and this contradicts the fact 
that V is bounded from above. Hence it must be H — > (faster than as t — > +00. 

Since if = K(b 2 /2, then (faster than i^ 1 / 2 ) as well when t — > +00. The 

asymptotic state for a potential bounded from above is one with H — > and 0^0. 
Eq. (|2.14j) gives the asymptotic relation i?Q = kV/3 and the Klein-Gordon equation 
reduces, in this limit, to dV/d(j)\^ Q = and the asymptotic state is the de Sitter fixed 
point (H ,(f> ). The function 

L(H,<P) = (H-H ) 2 (3.2) 

defined on the entire phase space is a Ljapunov function - in fact L(H, <p) > V (H, </>) ^ 
(H ,<j)o), L(H ,(j) ) = 0, and dL/dt = n<p 2 (H — H ) < V(if,0) ^ (H Q ,<p ) because 
H tends to H while non-decreasing, hence H < H Vi. The system is stable and the 
attraction basin of the de Sitter attractor is the entire phase space. This excludes the 
possibility of a Big Rip or a sudden future singularity. 

If the potential has maxima, minima, or inflexion points, there will be fixed points 
of the system (J2.3|) - (j2.5|) . The character of these points determines the stability of the 
fixed point. One can assess the stability by considering the perturbed equations ([2.4)1 - 
()2.5|) for homogeneous perturbations that depend only on time, as done in j2Hl 139] • 
It is more meaningful to study stability against inhomogeneous perturbations which 
depend on both space and time. Then the complication arises that a gauge-independent 
formalism is required. The study of inhomogeneous perturbations is carried out in [3U] 
for generalized gravity theories described by the action 



S — Id x 



\f{<j>,R)-^-g ah V a <jN b <j>-V{<j>) 



(3.3) 



The result of 00] is that there is stability if and only if 

— VwWo ) + MM 



UJQ (l + 



2wo/fl(</>o) 



< . (3.4) 



Phantom cosmology is the special case /(</>, R) = R, u = —1, /^r = — 0, /r = 1 
and stability corresponds to V " < 0, where a prime denotes differentiation with respect 
to cf). The perturbation analysis of jlU] establishes that the de Sitter fixed point is stable 
against inhomogeneous linear perturbations if V(4>) has a maximum there. Here we 
recover this result in a different approach that does not rely on a linear perturbation 
analysis and it is valid to any order, extending and complementing the results of jlU] 
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applied to phantom cosmology. Moreover, [JU] does not draw conclusions about the size 
of the attraction basin, while here we establish that this attraction basin is the entire 
phase space. This agrees with, and extends, the result of Guo et al. |H] that, if a late 
time attractor exists in the phase space of phantom cosmology with bounded potential, 
it is unique. 

The phantom field settling in a maximum and producing a late time de Sitter attrac- 
tor is consistent with the literature selecting specific potentials V(<j)) and with general 
conjectures jHHHSlESl- Let us see some examples. Singh, Sami and Dadhich [33] con- 
sider the bell-shaped potential V((p) — Vq/ cosh (nacfi) attaining its global maximum Vq 
at = 0, and they find a late time de Sitter attractor. The same type of attractor is 
found by Carroll, Hoffman and Trodden ^3] for the Gaussian potential V((p) = Vq 1° ' . 

Hao and Li [39j consider the inverse power- law potential V(<f)) = Vq — o (0/0o) 2 f° r 

|0| — ( f } o\f^'> with maximum Vq at = 0. They find again a de Sitter attractor empty 
of scalar field (0 = 0) and with a truly constant cosmological constant A = Vq. 

Another situation in which the potential is bounded from above occurs when V(<f>) has 
a horizontal asymptote Vo approached from below as — > ±oo. Under these conditions 
the scalar field rolls up the slope of the potential but its (negative) kinetic energy is 
dissipated by friction - described by the term 3H<f) in eq. ()2.5|) - while the force term 
—dV/d<fi tends to zero approaching the asymptote, with 3i?0 ~ dV/dcj) and 0^0. 
As a result the motion stops (0 — > 0) and a de Sitter regime with constant and 
H is approached. An example is the potential studied by Sami and Toporensky [31] 



V{4>) = Vq 



1 - e~ c ^ 



which has the shape of an inverted bell with Vq as horizontal 
lors find a late time de Sitter attractor with cosmological constant 



asymptote. These aut 
A = V |3J. 

Guo et al. jH] consider the bounded potential V{4>) = Vq [1 + cos (0//)] and find a 
slow-climb solution in which settles in the maximum of the potential at = in a de 

Sitter regime with (if o ,0o) = ^a/ k ^o/3, 



4 Unbounded potential 

If the phantom potential is not bounded from above it would seem that the asymptotic 
state could be either a de Sitter regime, or a very different state, possibly a Big Rip, 
depending on the shape of V((f)) and its slope. Physically, the motion of 0(£) is analogous 
to the motion of a ball with negative kinetic energy falling up a hill under the force dV j d<p, 
damped by the term 3if0. If this friction becomes negligible in comparison with the 
force and inertial terms, 3H(p << dV/d(p,<j), one has a "slow-climb" regime analogous 
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to the slow-roll, potential-dominated, regime of inflation for a canonical scalar field 
[SH E]. This happens if the potential becomes sufficiently steep. If instead friction is 
comparable to the force term, motion could stop (H ~ 0, <fi ~ 0) corresponding to a de 
Sitter regime. If V is strictly increasing this situation is forbidden and the expansion 
is super-exponential: if V((p) is not bounded from above and is strictly increasing the 
solution of eqs. ^.^ -i TOj) cannot approach a de Sitter fixed point at late times, but H 
and (f) go to infinity. 

In fact, a de Sitter fixed point has V = and in order to approach a de Sitter fixed 
point it must be V — > 0, which is incompatible with our assumptions. Hence H and 4> 
go to infinity. 

This result explains the examples available in the literature, in which an unbounded 
potential with dVj d<p > never produces a late time de Sitter solution but rather gives 
a universe that expands faster than exponentially. Sami and Toporensky [33] and Guo 
et al. jH] consider the simple potential V(<p) = m 2 <f) 2 /2. Let us consider the region 
4> > where dV/d(j) is positive - there <fi and H approach constant values at late times 

as t — > +oo, and H, oc t asymptotically. The motion becomes potential-dominated 
at large times and the ratio of kinetic to potential energy m m ~J^ 2 — > so that 

V » <p 2 /2 asymptotically. Both H and cf> diverge and the scale factor a(t) = aoe a * 2//2 
diverges faster than exponentially. However, it takes an infinite time for <p and H to 
reach infinity and there is no Big Rip. The effective equation of state parameter 

P (b 2 + m 2 6 2 , . 

w = —= — -> -1 (4.2) 

p (f) 2 — m 2 (f) 2 

as t — ► +oo. The asymptotic equation of state approximates P = —p, which is the exact 
equation of state for de Sitter space, but the solution is not de Sitter space nor does it 
approach it. Asymptotically, — > and the friction and force terms balance each other, 
3H(f) ~ dV/d(j). The appendix is devoted to this example and shows that (j4.1|) is an 
attractor. 

Sami and Toporensky jHH also consider power-law potentials V(<fi) = Vo4> a (with 
a > 0) and they find a "slow-climb" regime characterized by 
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if a < 4. For a = 4 they find exponential growth of both and H, while if the power-law 

2 

potential is steeper (a > 4), they find a Big Rip solution with <f>(t) ~ (to ~ t) 4 ~ a . This 
result is confirmed by Guo et al. [H] and agrees with the qualitative argument on the 
steepness of the potential V((j)). 

Hao and Li consider the exponential potential V(4>) = Voexp (c0) (with c > 0), 
which is unbounded and has dV/dcf) > everywhere. They find an attractor with 
equation of state parameter w < — 1 which makes the Big Rip unavoidable. In [35] the 
same authors consider again the exponential potential V(4>) = Voexp \\/ZkA(4> — 0o)] 
with Vo, A > finding a Big Rip attractor in parametric form (a(4>),t((f>)). By eliminating 
the parameter one finds the scale factor 



^ 2^+2 ) 

In simple terms, the effective equation of state parameter is constant, w = —(1 + A 2 ), 

2 

and the scale factor can be written as a oc (t* — t) 3 (»+ 1 ) . It is well known that a spatially 
flat FLRW universe filled with a fluid with constant equation of state P = wp has scale 

2 

factor a oc £3(i+«o 

It seems intuitive that a potential V(<p) unbounded and steeper than the exponential 
potential will always cause a sudden future singularity. This can be shown for extremely 
steep potentials: if V(<p) has a vertical asymptote the universe evolves to a future sin- 
gularity in a finite time. 

In fact, since V(0) — > +oo as (ft — > <fi c , it is dV/d(p > 0. V is unbounded and 
strictly increasing and, due to our previous result, cf) cannot tend to zero. The only two 
possibilities are that <p tends to a finite limit <fi c or that <p diverges. In both cases </> reaches 
the critical value <p c in a finite time and V and dV/dcf) diverge. <fi cannot have a finite 
limit C : in fact, if this happens, then H = K(p 2 /2 — > H c = n<f) 2 ./2 and (f) — > 0. Then the 
modified Klein-Gordon equation ()2.5|) yields 3H(p ps dV/d<p — ■> +oo asymptotically and 
if must diverge in a finite time, which is in contradiction with H ~const. Hence it must 
be 4> — > +oo and H = k0 2 /2 — > +oo. It may happen that both H and </> diverge or that 
they stay finite with their derivatives diverging at a sudden future singularity. Examples 
of this kind of singularities have been found in other models of dark energy based on the 
tachyonic field jH], the brane world [13], and cosmology with the Gauss-Bonnet term 

Out last result is consistent with the numerical examples of (121 El]; who find a 
sudden future singularity for potentials with a vertical asymptote. 



2vM 2 + 2 



+ t -t 



2 



(4.4) 
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5 Discussion and conclusions 



The dynamics of phantom cosmology are studied without assuming specific potentials 
V ((/)), but assuming that the phantom dark energy has already come to dominate the 
cosmic dynamics. In the literature there are also scenarios in which the late time state of 
the universe is dominated again by cold (dark and ordinary) matter with zero pressure 
and the phantom energy decays [39 : such scenarios are a priori excluded from our 
analysis. 

A clear and unified picture of the dynamics is obtained: the phase space is a two- 
dimensional surface in the ^H, 0, <pj space and is composed of two sheets. All the 

equilibrium points, which exist if and only if dV/dcf) has zeros, are de Sitter spaces located 
on the curve where the two sheets touch each other - they are given by eqs. ([2.110 . H 
is always increasing except when <fi vanishes, and there are no periodic orbits. In the 
presence of a bounded potential V(0) the universe has a late time de Sitter attractor 
whose attraction basin is the entire phase space. If V(<f>) is unbounded there cannot 
be a de Sitter attractor and the late time cosmic expansion is either super-exponential 
with infinite lifetime or else the universe ends in a Big Rip or other sudden future 
singularity. Unfortunately, in the general situation, we cannot provide a sharp criterion 
on the potential to discriminate between these two possibilities. Many results for specific 
potentials that are found in the literature are recovered and predictions are made for 
new scenarios featuring potentials with the shapes discussed. 

It is worth remembering that there are serious doubts on whether phantom fields 
can be stable [T31 [2H1 12HI and that the classical equations of motion considered have to 
be superseded by a semiclassical treatment near the Big Rip singularity. Indeed there is 
some evidence that semiclassical backreaction may avoid the Big Rip [49 . 

Finally, there are compelling arguments for a field with canonical kinetic energy to 
couple nonminimally to the Ricci curvature [5U[ ITT] . In particular, conformal coupling 
is associated with an infrared fixed point of the renormalization group [STJ and avoids 
causal pathologies [521- It is not clear whether the same arguments apply to a phantom 
field and force it to couple nonminimally (a nonminimally coupled phantom is considered 
in |4T[ |2*9~] ) . The slow-climb regime for a phantom would be altered by the inclusion 
of non-minimal coupling terms in eqs. ()2.3[) - (j2.5[) . similarly to what happens in slow- 
roll inflation [IxJj. The dynamics of a nonminimally coupled phantom will be studied 
elsewhere. 
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Appendix: massive phantom field 

Here we consider a phantom field with potential V{4>) = m 2 2 /2 as an example. The 
phase space is the two-dimensional surface ^H 2 +(f> 2 = m 2 (j) 2 . By using the dimensionless 
variables ^ ^ 

x = ^kH, y = —<p, 2 = — 0, (A.l) 

the equation of the surface can be written as y = ±\/x 2 + z 2 , which describes a cone 
with axis along the y-axis and vertex in the origin. The upper sheet is described by 
y = +\ / x 2 + z 2 and the lower sheet by y = — \/x 2 + z 2 ; the two sheets join on the plane 
z = along the lines y = ±x. Although we restrict our analysis to spatially flat (K = 0) 
universes, it can be shown that the cone separates trajectories belonging to an open 
(K = —1) universe lying inside the cone from orbits belonging to a closed (K = +1) 
universe lying outside the cone (cf. |K3j). 

The only equilibrium point satisfying eqs. (j2.11j) is the Minkowski space (x, y, z) = 
(0,0,0). Our results imply that H and go to infinity for any solution except the 
equilibrium point, and hence the equilibrium point is unstable. Since V"(0) = m 2 > 
the fixed point (which coincides with the global minimum of the potential) is unstable 
- see the discussion of Sec. III. 

This three-dimensional picture should be compared with fig. 3 of [H] showing a 
projection of the phase space (a two-dimensional cone) and of some orbits onto the 

(<ft, 0^ plane, corresponding to our (y, z) plane. A similar conical geometry of the phase 

space structure appears if one considers inflation realized by a massive scalar with the 
"right" sign of the kinetic energy density 

The trajectories of the solutions exhibit a late time attractor found analytically in 
[HI] (eq. (jHU). The authors of show numerically the convergence of many 

solutions to this asymptotic solution but do not prove analytically its attractor nature. 
The attractor can be obtained by requiring an asymptotic linear relationship between 
H and 0, as suggested by numerical integrations. By setting H = acfi one obtains, using 

2 



eq. (J2X3D, = 2a/ k. Eq. yields H w ±yj^amt, while eq. (g3D is identically 

satisfied and eq. (|2.5J) gives a = m^J k/6. Since H > the negative sign is rejected. The 
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asymptotic solution found in this way coincides with the one of 

(#*,</>*)= |yi,^*j • (A.2) 

The stability with respect to linear homogeneous perturbations is assessed by considering 
the perturbed Hubble parameter and scalar field H = H*+5H, (f> = <j)^+S<p, and inserting 
these expressions into the evolution equations ()2.3|) - (j2.5J) . It is convenient to consider 

the contrasts Ai = 5H/H* and A 2 = 50/0*. One finds A 1 = (d(f) - ^) and 

eq. (|2.5jl yields the evolution equation for A 2 

A 2 + - m + m 2 ?j A 2 + (m 2 - A 2 = (A.3) 

with asymptotic solution A 2 ^const. and Ai — A 2 (l — — jjy) ~ — constant. The 
solution (H*,<fi*) is stable and is an attractor, in agreement with the numerical results. 
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